ASYMPTOTIC AND EXACT PRICING OF OPTIONS ON 

VARIANCE 



MARTIN KELLER-RESSEL AND JOHANNES MUHLE-KARBE 



Abstract. We consider the pricing of derivatives written on the discretely 
sampled realized variance of an underlying security. In the literature, the re- 
alized variance is usually approximated by its continuous-time limit, the qua- 
dratic variation of the underlying log-price. Here, we characterize the small- 
time limits of options on both objects. We find that the difference between 
them strongly depends on whether or not the stock price process has jumps. 
Subsequently, we propose two new methods to evaluate the price of options on 
the discretely sampled realized variance. One of the methods is approxima- 
tive; it is based on correcting prices of options on quadratic variation by our 
asymptotic results. The other method is exact; it uses a novel randomization 
approach and applies Fourier-Laplace techniques. We compare the methods 
and illustrate our results by some numerical examples. 



1. Introduction 

Consider a discounted asset S = So exp(X) and a time-interval [0, T] subdivided 
into n intervals of equal length with boundary points tj = j— for j = 1, . . . , n. The 
corresponding (annualized) realized variance of X over [0, T] is then defined as 

n 1 n 

(1.1) RV n x (T) = - 5>g(,VVx) 2 = T E ~ X ^f- 

j=i o=l 

There exists a considerable number of financial instruments that are based on real- 
ized variance as an underlying (see, e.g., [3 [8] or [TTJ Chapter 11] for an overview). 
Well-known examples are variance and volatility swaps, as well as puts and calls on 
realized variance with payoffs (K — RVr) + resp. (RVt — K) + . By market conven- 
tion, the length tj — tj—± of a single interval typically corresponds to one business 
day for these derivatives (see [SUB])- For puts and calls, the notion At-The- Money 
(henceforth ATM) refers to choosing the strike K equal the to swap rate, which in 
turn equals the expectation E [RV* (T)] under the pricing measure. 

Given a stochastic model for S resp. X, the standard approach to pricing options 
on realized variance is to approximate realized variance by 

(1.2) RVX(T)nhx,X] T , 
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where [X, X] is the quadratic variation of the log-price X. This approximation 



is motivated by the fact that - for fixed T - realized variance (1.1) converges to 
h\X, X\i in probability as the number of subdivisions n tends to infinity (cf., e.g., 
[131 Theorem 1.4.47]). The advantage of this approach is that for many stochastic 
processes, the quadratic variation is a well-studied object. For example, as recently 
shown by |16j . the characteristic function of the quadratic variation in any affine 
stochastic volatility modeQcan be computed as the solution of a generalized Riccati 
differential equation, such that in many cases methods based on Fourier-Laplace 
inversion (cf. 6, 22J) can be applied to compute option prices efficiently. Moreover, 
using quadratic variation in place of realized variance, allows - at least in diffusion 
models - for elegant replication arguments, such as the the representation of a 
variance swap as an infinite portfolio of European options (see |20j). 



The quality of the approximation (1.2), or more precisely, the speed of conver- 



gence of realized variance to quadratic variation as the number of subdivisions n 
increases has been studied extensively in the econometric literature. Barndorff- 
Nielsen and Shephard [2J, for example, obtain a central limit law for the difference 
between realized variance and quadratic variation, scaled by the square root of n, 
which holds in a large class of stochastic volatility models (compare also [14] for 
more general results). 

However, an important difference from the econometric setting to the setting of 
variance options is that the sampling frequency n for such options can not be chosen 
freely, but is determined by market convention. As mentioned above, daily sampling 
is the typical case. From an econometric point of view, this sampling frequency 
would most likely be considered insufficiently low to produce an acceptable estimate 
of quadratic variation over most reasonable time horizons T. For variance options, 
though, not pointwise estimation under the physical measure is the goal, but the 
accurate risk-neutral pricing and hedging of options with specific payoffs. Therefore 
it is not clear to what extent findings from econometrics can be transferred. For 
this very reason several articles have considered the quality of the approximation 



(1.2) purely from the point of view of option pricing. Biihler [5 and Sepp 



resp. Broadie and Jain [3] find via Monte-Carlo simulation resp. analytically that 



the approximation (1.2) for daily sampled realized variance works very well for 
claims with linear payoffs, like variance swaps. On the other hand, Biihler [5] 
observes that "while the approximation of realized variance via quadratic variation 
works very well for variance swaps, it is not sufficient for non-linear payoffs with 
short maturities. The effect is common to all variance curve models ( or stochastic 
volatility models, for that matter).^ In particular, he presents some numerical 
examples based on call options on realized variance in the Heston model, that 
indicate that the approximation by quadratic variation notably diverges from the 
true value for short maturities (cf. p. 128]). This leads to the following questions 
considered in the present study: 

(1) To what extent is it indeed true that quadratic variation is not a good 
proxy for realized variance, when pricing short-dated options with non- 
linear payoffs? 

(2) How can options on the realized variance be valuated more accurately? 
The remainder of the article is organized as follows. We first consider exponential 

Levy models in Section[2]and compare the small-time limits of options on quadratic 
variation and on realized variance. Based on the results we propose a new method 
to approximatively evaluate prices of options on realized variance. In Section [3] 
we show that the results on exponential Levy models can be transferred without 



lr The class of affine stochastic volatility models includes exponential Levy models, the Heston 
model with and without jumps, and many stochastic time-change models. 
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modification to general semimartingale models. In Section[4]we propose a second - 
exact - pricing method for options on realized variance, which is based on Fourier- 
Laplace methods. We present numerical examples and compare the two methods 
in Section [5] and then conclude with some suggestions for future research. 

2. Small-time asymptotics in exponential Levy models 

In this section, we derive the small-time asymptotics for options on variance in 
exponential Levy models. That is, we suppose the asset price process is modeled 
as S = Sq exp(X) for a Levy process X. The latter is assumed to be square- 
integrable, such that the variance swap rate is always finite. The Levy process X 
will be characterized through its Levy-Khintchine triplet (6, a 2 , F(dx)) with respect 
to the truncation function h{x) — x or, equivalently, by its Levy exponent, i.e., the 
function 

ip(u) = ub+ -u 2 a 2 + / (e ux - 1 - ux)F(dx), u E iR, 



2 

for which E [e uXt ] = exp(tip(u)). We can decompose X as 

X t = U + aWt + Lt, 

where W is a standard Brownian motion and L is an independent centered pure- 
jump Levy process. 

2.1. Options on quadratic variation. First, we consider the simpler case of 
options written on (annualized) quadratic variation. 

Theorem 2.1. Let X be a square- integrable Levy process with Levy-Khintchine 
triplet (b,o~ ,F{dx)) and suppose the payoff functions gx : R+ — > R, T > are 
continuous, uniformly bounded, and satisfy \ \gx — <?o| loo ~ ^ as T — > 0. Then 

HmE[ 3T (I[X,l] T )] = .go(^ 2 ). 

T— >0 

The proof is based on the following auxiliary result, which is shown in the ap- 
pendix. 

Lemma 2.2. Denote by L the pure-jump component of X . Then we have 

^[L,L] T -> a.s. as T -> 0. 

Proof of Theorem \2.1\ Evidently, 

\E[g T (±[X,X] T )] - 90 (<J 2 )\ 

< E [\g T {^[X,X] T ) -g (±[X,X] T )\] +E[\g (UX,X} T ) ^ g (a 2 )\] . 

By dominated convergence, the first term converges to zero, because we have \ \g? — 
5o||oo — > as T — > 0. Likewise, dominated convergence and the continuity of 50 
imply that the second term also converges to zero as T — > 0, because ^[X, X]t — 

□ 



rJ - 1 - }^[L,L]t — > o~ 2 a.s., by Lemma 



2.2 



2.2. Options on Realized Variance. The analogue of Theorem |2.1| for options 
on the discrete realized variance (1.1) reads as follows: 

Theorem 2.3. Let X be a square-integrable Levy process with Levy-Khintchine 
triplet (b,a 2 ,F(dx)) and suppose that the payoff functions g n ,T '■ — > K, T > 
0,rt G N are uniformly bounded and satisfy \\g n ,T ~ ffn,o||oo ~^ as T — > for each 
n. Then 

(2.1) hm E [gnAEV* (T)] = E [g n , Q (Y n )} , 

where Y n has gamma distribution with shape parameter n/2 and scale parameter 
2a 2 /n. 
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Note that the distribution of the limiting random variable Y n is determined 
solely by the number n of sampling dates and the (annualized) variance a 2 of the 
continuous part of X. The jump part of the process can only affect the value 
indirectly through the payoff functions g n ,Ti which usually depend on the swap 
rate and hence on the specification of X. For the most important case of calls 



and puts on variance, it will turn out in Section 2.3 below that the final result 
indeed depends also on the jump part, namely through the (annualized) variance 
V = J x 2 F{dx) of the jump measure. For ATM calls and puts, the small-time limit 
will be determined completely by the value of v 2 . 

To prove Theorem|2.3| we need the following two lemmas. The proof of the first 



one is similar to the proof of Lemma 2.2 and can be found in the appendix 



Lemma 2.4. Denote by L the pure-jump component of X. Then we have 

RV„(T) -> a.s. as T ->• 0. 

Lemma 2.5. Let C = o~(L t ;0 < t < T) be the a-algebra generated by the pure- 
jump component L of X . Then conditionally on C, the rescaled realized variance 
^RV* (T) follows a non-central chi-square distribution, with n degrees of freedom 
and noncentrality parameter 

A(T)=W 2 T+5l r+ 4<(T). 

Proof. Conditionally on £, the independent random variables X tj — X tj _ 1 , j = 
1, . . . , n are normally distributed with 

E [X t . - X t ._ t | £] = IT In + L t] - L t ._ 1 , 
Var [X tj -Xt-^C] = a 2 T/n. 

The rescaled realized variance ^RV*{T) is therefore non-central chi-square dis- 
tributed with n degrees of freedom, conditionally on C The noncentrality param- 
eter of the distribution is given by 



j^N & v[X tj -X t ^\C] 



b\ 2 _ 26 



- ^+^ + -^£(L,-L^, 

7 3=1 



as claimed. □ 



Proof of Theorem \2.3\ Denote by 

O-n/2 

(2-2) f n {x) = ———x n/2 ~ 1 e~ x/2 , x>0, 

T(n/2) 

the density of the central chi-square distribution with n degrees of freedom. The 
density f n ,\{x) of the non-central chi-squared distribution with n degrees of freedom 
and non-centrality parameter A > can be expressed as an infinite weighted sum 
of densities of central chi-square distributions: 

00 -A/2A 

fnA X )=J2 T^/n+2i(*), X>0. 

i=0 
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Finally, note that -/„ (-), with c > is the density of a gamma distribution with 
shape parameter n/2 and scale parameter 2c. Using Lemma 2.5 we have 

E [g n ,T(RVX (T))] = E [E [g n>T (RV* (T)) \ C]) 
t2 x / e -A(T)/ 2(A(T)/2) i 



(2.3) 



E 



, OO 

E 

i=0 



9n,T 



-fn+2i(x) dx 



Let M be the counting measure that assigns mass 1 to each integer in No- The last 



term in (2.3) can then be regarded as integrating the function 

' e -A(T)/2 (A(T)/2 y 



h ntT (u,x,i) = g n , T 



fn+2i(x) 



with respect to the product measure P <g> da; ® M. We want to evaluate this integral 
as T — > 0. By Lemma 2.4 A(T) — > P-almost surely. Hence we may assume 
that T is small enough to ensure A(T) < 1, almost everywhere with respect to 
¥(&dx®M. In this case we can estimate the integrand using the explicit form (2.2) 
of the chi-square density: 

(x/4Y 



\hn }T (u,x,i)\ < \\g n ,q 



Or/2) 



i/2-l x / 2 



Summing the right-hand side with respect to the counting measure M we get 
||<?n.,:r|| (x /2) n / 2 ~ 1 e^ x / i , which is P <S> cfe-integrable. Therefore, dominated con- 
vergence allows us to interchange limit and integration, and we obtain 



limE[ Sll!T «(T))] =E 



Ehm h n r(w, x, i)dx 

i=0 



'n,0 ( —X I fn(x)dx 



Using the fact that a scaled chi-square distribution is a gamma distribution, the 
result follows. □ 

Having determined the small-time limit of option prices on both quadratic vari- 
ation and realized variance, we now consider the difference between the two. 

Definition 2.6. Assume that the limiting payoff g is the same for all options, i.e., 
there exists a function g : R + — > K such that \\gr — gW^ — > and \\g n ,T — gW^ — > 
as T — > for all n G N. Then we define 



lim (E [g n!T (RV*(T))] - E [gr$[X, X] T )] ) , 



(2-4) A n (g) 
and call A n (g) the discretization gap. 

Corollary 2.7. Suppose the prerequisites of Theorems \2.1\ and are satisfied. 



Then the discretization gap A„(g) from Definition 2.6 is given by 
(2-5) A n (g) = E [g(y n ) - g(a 2 )] , 



where Y n is a gamma- distributed random variable as in Theorem 2.3 If, in addition, 
the function g is convex, A n (g) has the following properties: 

(a) A n (g) > for all n£N, 

(b) A„(<?) — if and only if u 2 = or g is affine-linear, 

(c) n i— > A n (g) is decreasing in n and converges to as n — > oo. 

Proof. Since g is convex and E [Y n ] — a 1 , Jensen's inequality yields E [g(K„)] > 
E [<7(cr 2 )] and Q follows. Equality clearly holds if a 1 = or g is affine-linear, 
which yields the 'if '-part of Q. For the 'only if '-part assume that a 1 > and that 
g is not affine-linear. Then g is strictly convex at least on some interval (a, b). Since 
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a 2 > the interval (a, b) has strictly positive measure under the law of Y n and the 
strict Jensen inequality implies that A n (g) > 0, completing the proof of Q. By 
[281 Example 1.5.1e] the gamma distributed random variables Y n are decreasing in 
the convex stochastic order. In particular, n 4 E [sC^n) — g( (j2 )] is decreasing, 
too. Finally Y n converges to a 2 in distribution as n — > oo by elementary properties 
of the gamma distribution; hence linin^oo A n (g) — 0, showing Q. □ 

The above corollary has some very interesting implications. 

• Assertion ([a]) shows that - at least asymptotically for small maturity - an 
option on quadratic variation is always cheaper than the option on realized 
variance with the same payoff, given that the payoff is convex. 

• Assertion Q shows that the difference between the two option prices - the 
discretization gap - vanishes in two cases: The first case is when the payoff 
is linear; this confirms the observation of Biihler quoted in the introduction, 
and explains why for variance swaps realized variance can be substituted 
by quadratic variation even for short maturities. The second case in which 
the discretization gap vanishes is for Levy processes without a diffusion 
component, (i.e., with a 2 = 0). This suggests that in a pure-jump Levy 
model, quadratic variation should be a good proxy for realized variance, 
even when pricing short-dated options with non-linear (convex) payoffs. 
This assertion is confirmed by our numerical examples in Section [5] 

Assertion Q, finally, is also quite intuitive. Since realized variance converges to 
quadratic variation as n — > oo, also the discretization gap should vanish in the 
limit. It does, and in fact it does so monotonically in n. 

2.3. Applications to put and call options. Let us now examine the important 
special cases of put and call options. More specifically, denote by Vt the variance 
swap rate, and consider the payoffs 

x H> (kVr — x) + resp. x h-> (x — kVr) + 

for puts resp. calls with relative strike value k > 0. Setting k = 1 yields ATM 
options. If the realized variance is approximated by quadratic variation, the swap 
rate is given by 

V T := E [±[X,X] T ] = a 2 + v 2 , 



where v 2 = J x 2 F(dx). We may apply Theorem 2.1 to the payoff 

g T {x) = g (x) = (£E [[A, X] T ] ~ x) + = (k(a 2 + v 2 ) - x) + , 

to obtain the small-time limit for put options. For call options we use put-call 
parity: The difference between a call and a put option with relative strike k is 
(1 — k) times the swap rate. 

Corollary 2.8. Let X be a square-integrable Levy process with Levy-Khintchine 
triplet (b,a 2 ,F(dx)), and define v 2 = J x 2 F(dx). Then the following holds: 

(a) Put options on quadratic variation satisfy 

(2.6) YimE[(kV T - ±[X,X] T )+] = [a 2 {k - 1) + v 2 k] + , 

(b) Call options on quadratic variation satisfy 

(2.7) hm E [L [X, X] T - kV T ) + ] = v 2 + [a 2 (1 - k) - v 2 k] + . 

Note that in the ATM case both limits coincide and are equal to v 2 . In other 
words, the small-time limit of the price of an ATM option on quadratic variation 
is equal to the second moment of the jump measure and vanishes precisely in the 
Black-Scholes model. In the out-of-the-money case (puts with k < 1 and calls with 



ASYMPTOTIC AND EXACT PRICING OF OPTIONS ON VARIANCE 



7 



k > 1) the limit may be zero even when jumps are present. 

When the discretely sampled realized variance is used, calculations are a bit 
more involved. In this case, the swap rate is given by 

(2.8) V3? := E [RV r f(T)] = (cr 2 + v 2 ) + b 2 T/n. 

In particular V$ := UniT-»o Vt = (cr 2 + v 2 ). For put options Theorem 2.3 can be 
applied directly for call options we use again put-call parity. 

Corollary 2.9. Let X be a square-integrable Levy process with Levy-Khintchine 
triplet (6, a 2 , F(dx)), and define v 2 = J x 2 F(dx). Then the following holds: 

(a) Put options on realized variance satisfy 

2 \ / 2 



(2.9) UmE [{W£ - RV n (T))+] = <J 2 Q k , n (~j + (a 2 (k-l)+v 2 k)R Kn 

(b) Call options on realized variance satisfy 

(2.10) limE [(RV n (T) - kVp) + ] = 

= v 2 + a 2 Q Kn (^j + (a 2 {k - 1) + v 2 k) S.R k>n (j^j 

The functions Qk,n( r ) resp. Rk,n(r) are strictly decreasing resp. increasing func- 
tions on [0, oo ) given by 

n M 2/n (n k{\ + r) \ n/2 7 (n/2, k(l + r)n/2) 



r(n/2) \2 exp(fc(l +r))J ' , w T(n/2) 

whereT{x) denotes the (complete) gamma function and j(n, x) the lower incomplete 
gamma function. 



Note that it follows from Corollary 2.7 that both (2.9) and (2.10) are decreasing 
functions of n, which is illustrated by the numerical examples in Section [5] Also 
note that contrary to options on quadratic variation, the limiting value of a put 
or call option on realized variance is never zero, apart from the trivial case of 
deterministic X. 

Having derived the small-time limit for prices of puts and calls, both on realized 
variance and on quadratic variation, we can now consider the difference between 
the two, i.e., the discretization gap introduced in Definition |2.6| Simplifying the 
notation a bit, we write 

(2.11a) AP fc ,„ = Hm (E [(kV? - RV*(T))+] - E [(kV T - ±[X,X\ T ) + ]) 
(2.11b) AC\„ = Hm (E [(RV*(T) - kV£)+] - E [(±[X,X\ T - kV T )+]) , 

which are the discretization gaps for put and call payoffs respectively. The following 
results can be derived from Corollary 2.7 or simply by combining Corollaries 2.. 
and 12^91 above. 



Corollary 2.10. Let X be a square-integrable Levy process with Levy triplet (6, cr 2 , F(dx)). 
Set v 2 = J x 2 F{dx) and define the functions Qk,n(r) as well as Rk, n {r) as in Corol- 
lary^^ Then the following holds: 

(a) For a 2 (k — 1) + kv 2 > 0, we have 

AP fei „ = ACV = a 2 Q Kn ( V -^\ + (a 2 (k - 1) + v 2 k) [Ru n ( ?L) - l) . 



8 



MARTIN KELLER- RESSEL AND JOHANNES MUHLE-KARBE 



(b) For a 2 (k — 1) + kv 2 < 0, we have 



APfc : „ — ACk,n — 0- 2 Qk,n 



- 1) 



k)R 



From Corollary 2.7 it follows that AP n ^ and AC n .k are always positive and 
vanish if a 2 = 0, i.e., in a pure-jump model. Suppose now that, in a certain 



Levy model, we can easily calculate the prices of put and call options on quadratic 
variatiorj^J Then the expressions for AP n ^ and AC n ^k can be used as correction 
terms to obtain an improved approximation for the price of the corresponding option 
on realized variance: 



Approximation 2.11. Let AP n ^ and AC n ^k be given by Corollary 2.7 Then the 



price of a put resp. call on quadratic variation can be approximated by 
(2.12a) E [(kVi 1 - RV n {T)) + ^ « E \{kV T - i[A, A] T ) + ] + AP n , k , 

(2.12b) E [(RV n {T) - fcV?) + ] ^E[(±[X,X} T -kV T ) + ] + AC n , k . 

Remark 2.12. An approximation of similar type has been proposed by Sepp |27j 
for the Heston model. 

These approximations are exact in the limit T — ► (by definition of A.P n ^ 
and AC Ut k), and in the limit n — > oo (since the A-terms vanish, and realized 
variance converges to quadratic variation). We can therefore expect (2.12| to be 



good approximations for all maturities. The numerical results in Section [5j confirm 
convincingly that this is the case. 

Similar approximations can of course be constructed for general payoffs g, using 



the correction term A n (g) from ( 2.5 1. According to Corollary 2.7 the discretization 
gap A n (g) vanishes for linear payoffs, such that it can be interpreted as a convexity 
correction that corrects the basic approximation (1.2 1 depending on the convexity 
of the payoff. 



So far, we have confined ourselves to Levy models. It is a natural next step 
to examine whether our findings remain true when passing to more general asset 
price models incorporating, e.g., stochastic volatility. This is done in the following 
section. 



3. Small-time ASYMPTOTICS IN SEMIMARTINGALE MODELS 

In this section, we show that - under very mild conditions - the small-time 
asymptotics of options on variance for general semimartingales coincide with those 
of a suitable Levy approximation. Thus, the results derived in the previous section 
can be transferred directly and no new phenomena arise. Throughout, we suppose 
that the log-price is given by 

(3.1) dX t =b t dt + o- t dWt + K t {x)*(N(dt,dx)-F{dx)dt), A = 0, 

for a standard Brownian motion W, a Poisson random measure N(dt,dx) with 
absolutely continuous compensator F(dx)dt (cf. [T51 Section II. 1] for more details), 
and predictable integrands b, a, k. To ensure that the log-price process A is a 
well-defined square-integrable semimartingale, we assume that 

(3.2) £ (e[&?]+E[o?]+E 



K 2 (x)F(dx) 



dt < oo. 



'We discuss in Section 4 how - and in which models - this can be done. 
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Remark 3.1. The processes of the form (3.1) comprise essentially all semimartin- 
gales with absolutely continuous characteristics, i.e., without fixed times of discon- 
tinuity (see [T^l Theorem 14.68(a)] for more details). 



In any reasonable application, (|3.2| will imply 
(3.3) 



E 



K^(x)F(dx) 



< oo. 



In this case, the process X can be approximated for small t by the square-integrable 
Levy process 

(3.4) dX t = b dt + aodW t + Ko(x)*(N(dt,dx)-F(dx)dt), X = 0, 

obtained from X by "freezing" the coefficients of X at time zero. Subject to weak 
regularity assumptions on the coefficients b, a, and k, we then have the following 
small-time approximation results closely related to O Proposition 2.1]. 

Lemma 3.2. Suppose ( |3~2"| and ( pT3| hold. Then if E [|of — Oq|] -> and 
E [/ \k 2 (x) - K%(x)\F(dx)] -> fort^r 0, we have 

\imE[\±[X,X] T - ^[X,X} T \] =0. 

Proof. By [HI Theorems 1.4.52 and II. 1.8] and Fubini's theorem, we have 
E[\[X,X] t -[X,X T \] 



< f \[\a 2 -a 2 \]dt + f 
Jo Jo 



E 



\ K 2 (x) - K 2 (x)\F(dx) 



dt. 



Hence the assertion follows from the regularity assumptions on a and n. 



□ 



Lemma 3.3. Suppose ( |3.2[ ), ( |3.3[ ) hold and assume thatE [(6 t — 0o) 2 ] — ^ as weH 
as E [(er t - cr ) 2 ] ^ and E [J (^(x) - K (x)) 2 F(dx)] -> /or f -> 0. T/ien, /or 
any neN, we Ziaue 



lim E 

T-i-0 



i?y n x (T)-i?v; A (T) 



v , 



= 0, 



for the realized variances RV* (T) of the log-price X and RV*(T) of the Levy 
approximation X. 

Proof. First notice that the inequalities of Cauchy-Schwarz and Minkowski imply 



(3.5) E 



3=1 



3=1 



< 



£(e [(X^-X^+EiiX^-X^f)] 



2MV2 



J = l 



x E 



1/2 



E 



* 2 

tj_i 



1/2 



E 



1/2 



■E 



X 



1/2 



For t < T, the Burkholder-Davis-Gundy inequality as in [3TJ Theorem IV. 
well as QH Theorems 1.4.52 and II. 1.8] yield 

l-T 



E[X 2 ] <C^ (E[b 2 ] +E[a 2 ] +E 



K 2 (x)F(dx) 



dt. 



for a constant C which does not depend on T. Arguing analogously for X instead 
of X, we obtain 



(3.6) E[X 2 f 2 - 



-E 



X 2 t 



1/2 



2 1 1/2 



-E 



x 2 t 



1/2 



0(T 1/2 ) asT^O. 
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Now notice that another application of the Burkholder-Davis-Gundy inequality and 
[321 Theorems 1.4.52 and II.1.8] shows that, for t < T, 

E [(X t - X t ) 2 } 



< c 



E 



(«t(x) - kq[x)) F(dx) 



dt. 



E[{b t -b ) 2 ] + E [((Tf — Co) 2 ] 

for a constant C independent of T. Under the stated assumptions, it follows that 
E [(X u - X u ) 2 ] 1/2 + E [(X u _, - X ti _J 2 )] 1/2 = oiT 1 ' 2 ) as T -> 0. 



Combined with (3.5) and (3.6|, this proves the assertion 



□ 



For Lipschitz continuous payoffs, it now is an immediate consequence of Lemmas 



3.2 and 3.3 that the small-time asymptotics for the semimartingale (3.1| and its 



Levy approximation fl3.4| coincide. Therefore they are determined by the formulas 



in Theorems 2.1 and|2.3[ respectively. 



Theorem 3.4. Let X be a semimartingale of the form (3.1) satisfying the prereq- 
uisites of Lemma \3.S\ resp. Lemma \3.3\ 

(a) Suppose that the payoff functions gx ■ — > T > are continuous, uniformly 
bounded, and satisfy \\gx — ffolloo — ► as T — > 0. Moreover, suppose that go is 
Lipschitz continuous. Then 

KmE[g T (±[X,X} T )] = g (a 2 ). 

(b) Suppose that the payoff functions g Ut T '■ — > M, T > 0, n S N are uniformly 
bounded and satisfy \\g n ,T ~ ffn.olloo ^0 as T — > for each n G N. Moreover, 
suppose that the g n .o are Lipschitz continuous. Then 

lim E [g n!T (RVX(T))] = g n , (Y n ), 

where Y n has gamma distribution with shape parameter n/2 and scale parameter 
2a 2 /n. 



Proof. We decompose 

E[\g T (±[X,X] T ) 



5o(cr 2 )| 



< E [\g T (±[X,X] T ) -go(i[X,X\ T )\] 
+ E[\g (±[X,X} T )-g (t[X,X] T )\ 
+ E[\g (±lX,X} T )~g Q (a 2 )\] , 



where X is the approximating Levy process (3.4). The first term on the right-hand 
side can be bounded by \\gx — ffolloo and thus goes to zero as T —> 0. The second 
term can be bounded by CEjjifX, X]t — ^[X,X]t\\i where C is the Lipschitz 
constant of g^. By Lemma 3.2 this term goes to 0, too. The third term also 
converges to 0, by Theorem 2.1 and the claim H) follows. Assertion lb]) is shown 
in the same way, substituting^ [X, X) T by RV*\T) and using Lcmmal3.3 as well 
as Theorem O □ 



Consequently, Corollary |2.7| on the discretization gap, and Corollaries |2.8[ |2.9| 
and |2. 10] on put and call options on variance also hold in the present semimartin- 
gale setting if we substitute a 2 = a 2 and v 2 = J n^{x)F{dx). In particular, we 
find that the discussion following Corollary |2.7| can be completely transferred to 
the semimartingale setting, and that - again in an asymptotic sense - Biihler's 
statement quoted in the introduction holds generically for semimartingales with 
non-vanishing diffusion component, where Co > 0. More specifically, for a contin- 
uous semimartingale the small-time limit for ATM options on quadratic variation 
is zero, but its counterpart for realized variance is not. For semimartingales with 
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jumps and a continuous martingale part, both limits are non-zero, but there is a 
non-trivial discretization gap. In pure-jump models however, where a = 0, the two 
small-time limits coincide, suggesting that quadratic variation should be a good 
approximation even for short maturities. 

As an illustration, we show now how the prerequisites of Theorem |3.4| can be 
verified in some applications. For the sake of clarity, we do not strive for minimal 
conditions. 

Corollary 3.5 (Levy driven SDEs). Let f : M. — >• R be bounded and Lipschitz 
continuous and let L be a Levy process with Levy exponent ip L , whose Levy measure 
F L (dx) has bounded support. Then there exists a unique strong solution Y to the 
SDE 

dY t = f(Y t _)dL t , Y Q e K. 
Moreover, S — Soexp(X) is a martingale and the prerequisites of Theorem 3.4 are 
satisfied for the process 

dX t =dY t -Tp L (f(Y t _))dt. 

Proof. The first part of the assertion follows from the standard existence and 
uniqueness theorem for SDEs as in [5TJ Theorem V.6], because / is Lipschitz. 
Since the support of the Levy measure F L (dx) is bounded, L has finite moments 
and exponential moments of all order by [Ml Corollary 25.8 and Theorem 25.17]. 
Hence S is a local martingale by Ito's formula and the true martingale property is 
a consequence of [13 Proposition 1.4.50(c)]. Now notice that by definition of ip L , 



the process X is of the form (3.1 1 with 

b t = -^f(Y t .)- J(ef<*-> 



1- f(Y t ^)x)F L (dx), 



a t = Vc L f(Y t _), K t (x) = f(Y t _)x, 

for the Levy-Khintchine triplet (b L ,c L , F L (dx)) of L. Denote by M > the max- 
imum of the Lipschitz constant and the uniform bound for the function /. Then 
since Y t = Y t - a.s., for each t, we have 

E[\<rf-a 2 \] <2M 2 c i E[|r t -y |], 



E 



k?0) 



kI(x)\ F L {dx) 



< 2M Z 



x 2 F L (dx)E[\Y t -Y \}. 



Likewise, since / is bounded and F L (dx) has compact support, there exists a con- 
stant C > such that 



E[(b t -b Q ) 2 ] <2M 2 c L E[\Y t -Y \] + 2M 2 C J x 2 F L (dx)E [\Y t - Y \] 



Finally, 



E[(a t -a Q ) 2 ] <2M 2 c L E[\Y t -Y Q 



E 



K - Ko ) 2 F L (dx) 



< 2M 2 / x 2 F L {dx)E [\Y t -Y \}. 



Y is right-continuous. Combined with the Burkholder-Davis-Gundy inequality, 
similar arguments as above show that Y is also bounded in L 2 on any finite interval. 
Hence Y t — > Yq in L 1 and it follows that the conditions of Theorem 3^ are satisfied. 

□ 



The next corollary of Theorem |3 ,4| covers many stochastic volatility models from 
the empirical literature as, e.g., the ones of Bates [3], where a 2 follows a square-root 
process, and of Barndorff- Nielsen and Shephard 0, where a 2 is given by a Levy 
driven Ornstein-Uhlenbeck process. 
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Corollary 3.6 (Homogeneous jumps). Let n t (x) — Kq(x) be time-homogeneous, de- 
terministic, and such that J n^{x)F{dx) and J KQ ( X - )>1 e 2K °^ F(dx) are finite. Then 
for 

the stock price S = Soexp(X) is a local martingale and the conditions of Theorem 
3.4 hold, if a 2 is right- continuous and bounded in L 2 in some neighborhood of zero. 



Proof. Evidently, (3.2 1 is satisfied for sufficiently small T under the stated assump- 
tions. Moreover, Ito's formula shows that S is a local martingale. The regularity 
conditions on n in Thcorem |3.4| are trivially satisfied. The ones for a follow, because 
the processes a 2 and (a — a ) 2 are uniformly integrable and hence continuous in 
L l . ' □ 

4. Exact pricing methods for options on realized variance 

In the previous sections, we examined the small-time limits for options written 
on the quadratic variation and on the discretely sampled realized variance. We 
have also proposed a method, Approximation |2.11[ to approximate the price of 
an option on realized variance, given that the price of the corresponding option 
on quadratic variation is known. We first recall in this section how to compute 
these prices efficiently using Fourier-Laplace methods. We then propose a new ran- 
domization approach, that allows to use similar methods to directly determine the 
exact price of an option on realized variance in exponential Levy models, without 
the use Approximation [2 . 1 1| The two methods will then be compared numerically 
in Section [5j Throughout the section we assume that the log-price X follows a Levy 
process with the same properties as in Section [2] 

4.1. Option pricing using integral transform methods. We first recall how 
to price European-style options using the integral transform approach of |2"2"] . 
The key assumption is the existence of an integral representation of the option's 
payoff function / in the following sense: 

f(x) = / p(z)e- zx dz, 

J R—ioo 

for p : C — > C and R > such that v i-> p(R + iv) is integrable. 
Example 4.1. For a put option we have 

^ nR+ioo e Kz ^ roo / e K(R+iv) 



for x > and any R > (cf., e.g., Corollary 7.8]). 

In view of Fubini's theorem, the valuation of options which can be represented 
like this boils down to the computation of the Laplace transform of the underlying. 
E.g., for the put on quadratic variation we have 

1 r<>° / KT(R+iv) \ 

E[(K-±{X,X} T )+]=— Re f - E [exp {-{R + iv)[X, X] T )]j dv. 

Using the put-call parity (x — K) + = x — K + (K — x) + , this leads to the analogous 
formula 

E[(±[X,X\ T -K)+] 

= E[i[X,X] T ] -K+ — Re f E [exp (-(R + iv)[X, X] T )]\ dv 
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for calls on variance, provided that [X, X]t is integrable. Evidently, one just has 
to replace the normalized quadratic variation ^ [X, X] by RV X (T) to come up 
with the corresponding formulas for options on discretely sampled realized vari- 
ance. Summing up, it remains to compute the Laplace transforms of the objects of 
interest. 

4.2. Options on quadratic variation. For exponential Levy models, the qua- 
dratic variation process [X, X] also follows a Levy process (cf. [7] for the self- 
decomposable and Qj)] for the general case). More specifically, we have the follow- 
ing 

Lemma 4.2. Suppose the log-price X follows a Levy process with Levy-Khintchine 
triplet (b,a 2 ,F(dx)). Then [X, X] also is a Levy process and its Levy-Khintchine 
triplet is given by (a 2 , 0, F^ x,x ^ (dx)) relative to the truncation function h(x) = 0, 
where, 

p[x,x] (G) = J i G ( x 2) F ( dx ) } VG e B. 

Proof. See [H Lemma 4.1] . □ 

Combined with the Levy-Khintchine formula [Ml Theorem 8.1], this result im- 
mediately yields the required Laplace transform. 

Corollary 4.3. We have E [ e -^ x ' x ^] = exp[T^ x ^ (-u)}, for 

*™C->-(~ + /(^-W*>). M»)>o. 

Consequently, we obtain 

1 r°° / p KT(R+iv) , v \ 

(4,, E [<* - I[X,X W1 = ± I Re (^-pe— <-<-») * 

for puts on quadratic variation. Likewise, for calls on quadratic variation, 
(4.2) E[(i[X,X] T -if)+] 

(p \ i roo / KT(R-\-iv) \ 

S + / ^*)) - K + ^ I Re ») *, 

provided that J x 2 F(dx) < oo. Some examples from the literature where the Levy 
exponent t/)^ x,x ' of [X, X] can be computed in closed form are summarized in 
Section [5] below. 

4.3. Options on realized variance. In this section we develop a corresponding 



integral transform pricing method for claims on the discrete realized variance (1.1) 
in exponential Levy models, without any use of approximation. The method we 
present will in general be of similar computational complexity as the method based 
on quadratic variation, especially in cases where the Levy exponent ip of X is of 
more tractable form than the Levy exponent ipl x > x 1 of [X, X]. Such cases include 
for example subordination-based processes like the normal inverse Gaussian or gen- 



eralized hyperbolic process, see Remark 4.9 for more details 



As in the previous section, the crucial quantity is a Laplace transform, namely 



E 



exp ( -u^2(X t . -X tj _ x ) 



E 



exp(-uX|, /r 



where the equality is due to the independence and stationarity of the increments of 
the Levy process X. Consequently, if the Laplace transform of the squared process 
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is known, the price of, e.g., puts and calls on discrete realized variance can be 
recovered by an inverse Laplace transform as above. 

Our approach is based on the following identity: If Z is a normally distributed 
random variable, independent of X t , then using the characteristic function of the 
normal distribution it holds that 



(4.3) 



E 



-uXi 



E 



/2uX t Z 



E 



for all u G R+. Note that the first expectation is taken with respect to the law 
of the Levy process X t , the middle expectation with respect to the product law 
of X t and Z, and the final expectation with respect to the law of the normal 
random variable Z only. The exchange in the order of integration is justified by 
the Fubini theorem, and the fact that the integrands on the left and right hand 



side are bounded by 1 in absolute value. The benefit of formula (4.3 1 is to replace 
an integration with respect to the law of the Levy process - which is typically not 
known explicitly - by an integration with respect to a standard normal distribution. 
The characteristic exponent -0 which appears in the expectation on the right is in 
most cases analytically known and of considerably simpler form than the law of the 
Levy process. 



Let us remark here that the randomization approach of formula (4.3) can be 



extended to the Laplace transform of powers \X t \ p withp 6 (0, 2), and consequently 
to the discrete realized p-variation J2j l-^tj — -^j-il p °f a Levy process X. To 
this end, replace the standard normal variable Z by a symmetric a-stable random 



variable S p with parameters (a, 0, c,r) = (p, 0,1,0) (cf. 



Using that E 
(4.4) 

for all u G K+ 



a iwS p ~\ 



Theorem 14.15]). 



cxp(— \w\ p ) we obtain 



E 



-wlXtl" _ 



= E 



= E 



^(iSpii 1 '") 



Remarks on Laplace Inversion. The integral in formula (4.2) can be considered 



as inverting a Laplace transform by integration along a contour in the complex 
plane. There are many alternatives to this inversion method, see, e.g., |10j for an 
overview. Some of these methods only require knowledge of the Laplace transform 



on the positive real line, and thus seem tailor-made for formula (4.3) which holds 
unless further conditions are imposed - only on 
is probably the Post- Widder inversion formula 



The best-known such method 



(4.5) 



lim 



(-1)' 



1 



fi n X{n + l)/x), x>0, 



where / denotes the Laplace transform of a function /, and its n-th derivative. 
The Post- Widder method suffers from slow convergence and cancellation errors, 
and modifications such as the Gaver-Stehfest algorithm have been introduced to 
improve its performance. After implementing the Gaver-Stehfest algorithm and 
performing some numerical tests, we observed, however, that small errors in / — 
which invariably result from the evaluation of (4.3) - are strongly amplified by this 



method and lead to huge errors in /, probably due to the use of very high-order 



derivatives in (4.5). Moreover, as (TU] shows, inversion algorithms that evaluate the 
Laplace transform in the complex half-plane are in general numerically superior to 
algorithms that evaluate the Laplace transform only on the positive half-plane. For 



these reasons we decided to concentrate on the contour integration formula (4.2) 



and as a next step, to extend ( |4.3[ ) to the complex half plane Re (u) > 0. Note that 
this extension is also necessary for other more recent Laplace inversion algorithms 
like the methods proposed in pQ. 
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Extension to the complex half plane. Extending (4.3) to the complex half plane will 



not be possible without imposing some conditions on ip. The following is sufficient: 

Condition 4.4. The characteristic exponent ip has an analytic extension from the 
imaginary halfline jR+ to the sector 

3tt' 



A = |ti e C : - < arg(u) < 
Moreover, the extended function ip satisfies the growth bound 

(4.6) limsup Re W re ' 9 )) < forall G (f,f). 

r— >oo f 

Remark 4.5. An analytic function satisfying the growth bound ( |4.6[ ) is often said 
to be of order 2 and type in the sector A. 

An elementary symmetry argument shows that given the above condition, ip can 
also be analytically extended from the negative imaginary halfline — iR+ to the 
conjugate sector A. 

Lemma 4.6. Suppose ip satisfies Condition \4-4\ Then it can also be analytically 
extended from the halfline — iM. + to the conjugate sector A. Overall, ip has a unique 
extension to the hourglass shaped region A M = A U {0} U A, which is analytic on 
both A and A and satisfies the growth bound 

(4.7) limsup < for all 6 G (f , f ) . 

r— f ±oo r 

Proof. Suppose that ip satisfies Condition |4.4[ i.e., it is an analytic function defined 
on A. For u G A define ip(u) — ip(u). O n the imaginary axis, this definition agrees 
with the Levy-Khintchine representation of ip. The analyticity of ip on A follows 
directly, e.g., by verifying the Cauchy-Riemann differential equations. The growth 
bound on A M is an immediate consequence of the construction of the extension. □ 

We can now establish the central result of this section: 

Theorem 4.7. Let X t be a Levy process with characteristic exponent ip and let Z 
be an independent standard normal random variable. Then 



(4.8) E 



e -uX t 



E 



holds for all u on the positive real line. If X t satisfies Condition 4-4 then (4.8) 
holds for all u in the positive half-plane {u £ C : Re (u) > 0}, with ip denoting the 
unique analytic extension described in LemmaU 



Remark 4.8. The square root denotes the principal branch of the complex square 
root function with branch cut along the negative real line. 

Remark 4.9. For most Levy processes proposed in the literature, the Levy exponent 
ip can be computed in closed form. Hence, the evaluation of the Laplace transform 
of typically requires one numerical integration. The corresponding formula for 



the Laplace transform of [X, X]t in Corollary 4.3 is therefore simpler, if the integral 
J{e~ ux -l)F(dx)can be computed in closed form. However, even if this is possible 
as, e.g., for CGMY processes and the models of Merton and Kou, one usually has 
to employ special functions (cf. Section [5]) such that the numerical advantage is not 
too big. On the other hand, e.g., for NIG or generalized hyperbolic Levy processes, 
J [eT ux — l)F(dx) has to be evaluated using numerical quadrature, such that both 
formulas turn out to be of a similar complexity. 
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Proof. Let u € H+ := {iieC: Re (u) > 0}. The function u H> i\/2u (using the 
principal branch of the square root) is a single- valued analytic function on H+, 
mapping u to \/2juj exp (|(arg(w) + 7r)), and thus V.+ to A. With the normal 
random variable Z taking values in K it follows that iZ^/2u e A 154 . Let e > 0. Then 



(4.7) implies that there exists Mg > such that 
Retire* 6 )) < er 2 



(4.9) 
Thus 



Mg, for r G 



fc V4 ' 4 J 



exp 



2|u|e l(arsu+7r)/2 ))l < exp (e2\u\Z 2 + Mg) 



Note that Z 2 follows a chi-square distribution with one degree of freedom. The 
right hand side thus has a finite expectation of value (1 — 4e|u|)~ 1 / 2 e M '' , whenever 
|u| < l/(4e). Since e was arbitrary, it can be chosen small enough to satisfy this 
condition. We have shown that 



/(«) = E 



{|Z|<«} 



exists for all u e H + . Next we show that it is also analytic. Let Z n 
be a sequence of truncations of Z and define 

/ B («)=E[e*( <z -^" . 

Since "0 is continuous on A 1 * 1 , /„ — » / pointwise in Moreover, since the inte- 
grand is absolutely bounded for u in compacts, each f n is analytic in T-L+ (cf. [331 
Chapter 10, Exercise 15]). Let K be a compact subset of H+. On JC the bound 



(4.9) can be turned into a uniform bound 

Re(i/>(u)) < eR 2 + M, u € K, 

where R and M depend only on /C, and we again use the continuity of ip on A M . 
By the Cauchy-Schwarz inequality, we obtain 

\f{u)~ f n (u)\ 2 =E [exp(2Re(V( l Z^))) l { |z|>„}] 2 

< E [exp (e4i?Z 2 + 2M)] P(|Z| > n) 

= (1 - 8i?e)- 1 / 2 e 2M P(|Z| >n), 

for all it G /C. This shows that the convergence of /„ to / is uniform on compact 
subsets of T-L + . But analyticity is preserved by uniform convergence on compacts 
(cf. [221 Theorem 10.27]), such that / is analytic. We have now shown that both 
sides of (4.8) are well-defined analytic functions on T-L + . Since they coincide on the 
positive real line, they must coincide on all of ~H+, and the proof is complete. □ 



The following example shows that Condition |4.4| can not be reduced to analyticity 
in the sector A alone: 

Example 4.10. Let X t = N t — tj, where N t is a Poisson process with intensity 1, and 
7 = e — 1, such that e x is a martingale. The Levy exponent of this process is given 
by ip(u) = e u — 1 — wy. Clearly, tp has an analytic extension to the whole complex 
plane and, in particular, to the sector A. But Re (^(re 10 )) — e rcos ( 6 ' cos(r sin(6*)) — 



1 — r7COs(#), such that the growth condition (4.6) is not satisfied, e.g., in the 
direction 6 — 37r/8. Finally the formula (4.8) is not well-defined on the whole 
complex half-plane {u G C : Re (u) > 0}. Indeed, a tedious calculation shows that 
, +1 



Re E 



is infinite for, e.g., t = 1 and u = 3/8 — i/2, such that 



g e tip(iZV2u) 



does not exist. 
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Even though Theorem 4.7 fails in this simple case, Condition 4.4 holds for most 
Levy processes used in applications. 

Example 4.11. Condition |4.4| is satisfied for the following Levy processes. 

(1) Brownian motion: In this case, ip(u) = ^r(u 2 — u) is an entire function. 
Moreover, limsup,,^ Re{^(re t6 ))/r 2 = ^Re(e 2if> ) < for all 6 G (f , 2f ). 

(2) The Kou model: This jump-diffusion process corresponds to 



1 2 2 X + u A_ 
-a u 



u 



2 2/_|_ — u z/_ + u 

for A + , A_, > and /iel determined by the martingale condition 

ip{\) = 0. Again, ip obviously admits an analytic extension to A and, in 
addition, limsup^^ Re(^{re' e ))/r 2 = ^Re(e 2ie ) < for all 9 G (f , ^f). 
(3) The Merton model: For this jump-diffusion process, we have 



ip{u) = i-iu + — u 2 + A 



S 2 

exp I 7?i H u 2 I — 1 



fora>0, A,<5>0, 7GK and fj, G R determined by the martingale 
condition V-'(l) = 0. Consequently, V can be analytically extended to A. 
Furthermore, since Re(7re* e + ^-r 2 e 216 ') < for sufficiently large r, it follows 

that limsup^^ Re(ijj{re l8 ))/r 2 = ^-Re(e 2ie ) < for all 6 € (f , ^f). 
(4) NIG processes: In this pure jump specification, 



ip(u) = [iu + S(y / a 2 - f3 2 - v/a 2 - (fi + u) 2 ), 

where <5, a > 0, j3 € (—a, a) and fi is determined by the martingale condition 
ip(l) = 0. Once more, ip admits an analytic extension to A. Moreover, 
limsup,,^ Re(ip{re lf) ))/r 2 = for all 6 G (f , ^). 
(5) CGMF processes: These generalizations of the VG process correspond to 

ij}{u) = CT(-Y)((M - u) Y - M Y + (G + u) Y - G Y )), 

for parameters C, G, M > and Y < 2. In particular, -0 can be analytically 
extended to A and limsup,.^ Re(^(re ie ))/r 2 = for all 6 G (f , ^ ). 

Based on these examples and the above counterexample we conjecture that Con- 



dition 4.4 is related to the absolute continuity or smoothness of the Levy measure. 

5. Numerical Illustration 

We now consider three numerical examples. First, we take a look at the Black- 
Scholes model, then we turn to the pure-jump CGMY model. Finally, we also 
consider the jump-diffusion model of Kou. 

In the Black-Scholes model the distribution of realized variance is known explic- 
itly; it is the non-central chi-square distribution. Thus neither Approximation [2.11 



nor our exact method from Section [4] would be needed to compute prices of options 
on realized variance. Nevertheless, the Black-Scholes model can serve as a first 
test case to compare our two methods. Note that Approximation |2.11| amounts in 
this case to approximating the (true) noncentral chi-square distribution of realized 
variance by a central chi-square distribution. For a volatility parameter of a = 0.3, 
ATM call prices on realized variance, ATM call prices on quadratic variation and 
ATM call prices on quadratic variation corrected by the discretization gap ACi.. 
are depicted in Figure [T] for maturities up to 50 days. 

Evidently, the prices of ATM calls on realized variance cannot be distinguished 



from the convexity corrected approximation (2.12) by eye, showing that Approxi 



mation 2.11 performs remarkably well. 
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Figure 1. ATM call prices on realized variance, ATM call prices 
on quadratic variation resp. convexity corrected ATM call prices 
on quadratic variation in the Black-Scholes model. 



Also note that the prices of calls on the realized variance and of convexity cor- 
rected calls on quadratic variation converge to the prices of calls on quadratic 
variation (which are zero) for increasing maturity, but the rate appears to be even 
slower than in the results for the Heston model reported in [5]. In particular, using 
quadratic variation as a proxy for realized variance does not work well here, unless 



one uses the convexity correction (2.12) 



Next, we turn to the pure-jump CGMY process. By Section 4], 
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Y Y(l-Y) 
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where 
(5.1) 



I(k,v,t) 



2~ k — k/2 



T(k)U 



K 1 V 

2' 2' 47 



for the confluent hypergeometric U -function U. We use the calibrated (yearly) 
parameters 

C* = 0.3251, G = 3.7103, M = 18.4460, Y = 0.6029, 



from Table 1]. Corollary 4.3 Approximation 2.11 and Theorem 4.7 then lead 
to the results in Figure [2] 

Evidently, the three price curves cannot be distinguished by eye. In particular, 
quadratic variation seems to serve as an excellent proxy for realized variance at all 
maturities here, which drastically differs from the results reported for the Heston 
model in [5], and also from the results for the Black-Scholes model shown in Figure[lj 
This reflects the fact that the discretization gap vanishes for pure-jump models 
according to Corollary 2.10 In particular, the convexity correction term AC is 
zero in this case. 
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Figure 2. ATM call prices on realized variance, ATM call prices 
on quadratic variation resp. convexity corrected ATM call prices 
on quadratic variation in the CGMY model. 

As a third example we consider the model of Kou, which includes both jumps 
and a Brownian component. Similarly as above, the Levy exponent ^ X ' X ^ of [A, X] 
can again be expressed in terms of the confluent hypergeometric [/-function: 

tp [x ' x] (u) = <t 2 u + -u) - 1) - X-(iy-I(l,v-,-u) - 1), 



with the function / from (5.1). In this jump-diffusion model, we have a non- 



vanishing discretization gap by Corollary 2.10 Using the calibrated yearly param- 
eters 

a = 0.3, A + = 0.5955, v+ = 16.6667, A_ = 3.3745, v_ = 10, 
from |25l Section 7.3], we obtain the numerical results depicted in Figure [3j 
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Figure 3. ATM call prices on realized variance, ATM call prices 
on quadratic variation resp. convexity corrected ATM call prices 
on quadratic variation in the in the Kou model for a — 0.3. 



Again, first notice that the convexity corrected call prices on quadratic variation 
almost perfectly match the exact call prices on realized variance. Next, note that 
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whereas quadratic variation appears to serve as a much better approximation than 
in the Black-Scholes model here, the discretization gap is still significant for short 
maturities, unlike in the pure-jump CGMY model. By Corollary |2.10| we would 
expect the gap to shrink for smaller values of a. The effect of reducing a to 0.2, 
while keeping all other parameters the same is shown in Figure [4] and agrees with 
this prediction. 



0.06 




0.04 




+ 


RV 






o 


QV 


0.02 












X 


QV+AC! . 



Days 



Figure 4. ATM call prices on realized variance, ATM call prices 
on quadratic variation resp. convexity corrected ATM call prices 
on quadratic variation in the in the Kou model for a = 0.2. 



6. Conclusions and Outlook 

We have proposed two different methods to calculate prices of options on realized 
variance, that improve upon the standard approximation by quadratic variation. 
The first method, Approximation |2.11 was found to work very well for ATM 



options in three different Levy models with and without jumps. By the results in 
Section 3, it is also possible to apply the same approach to more general models, pro- 
vided that prices of options on quadratic variation can still be computed efficiently 
Hence one objective for future research will be to test its numerical performance for 
stochastic volatility models with and without jumps. For affine stochastic volatil- 
ity models (see for example [El H7J), a class which includes the Heston model, 
the SVJ and SVJJ models of [TT] and most time-change based stochastic volatility 
models, the results of [TOj could be used as a starting point. Since our approxima- 
tion method "freezes" the stochastic volatility at time zero, one would expect it to 
perform worse for stochastic volatility models. On the other hand, Sepp [37] has 
obtained encouraging results for the Heston model with a similar method. 

Our second method, the exact Fourier-Laplace approach from Section |4j could 
also be possibly extended to stochastic volatility models. Again the class of affine 
stochastic volatility models seems particularly suitable, since in such models the log- 
price X and the stochastic variance process V have a joint conditional characteristic 
function of the form 

(6.1) E[e uXt+wVt \T h ] =expU(t-h,u,w) + V h i>(t-h,u,w)+X h u\ 

for u,w £ iM.. It seems thus possible to use a conditional version of the identity 



(4.8) in each time-step between business days, and to use the special form of (6.1) 



to convert this conditional identity into a recursive algorithm for the computation 
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of the Laplace transform of realized variance. The delicate point is to find ana- 
lyticity conditions analogous to Condition |4.4| that allow to extend the identity to 
the positive half-plane {u € C : Re (u) > 0}. A rigorous analysis of the necessary 
technical conditions as well as an efficient numerical implementation for this case 
is also deferred to future research. 



Appendix A. Proof of Lemma 12.21 and 12.4 



Proof. We first show Lemma 2.2 Lemma 2.4 follows then after minor modifications 
of the proof. Let (T m ) m eN be a sequence converging to zero and let e > 0. For 
sufficiently small T m , we have 



(A.l) 



E 



[L,L]i 
T 



A 1 



< E 



[L,Lh 



A e 



T 

± rr 



Now notice that [L,L] is a Levy process with triplet (0,0,/ l.{x 2 )F(dx)) relative 
to the truncation function h(x) = 0, which can be used because [L,L] is of finite 
variation. Hence it follows from [13l Formula (5.8)] that 

~[L,L]T m Ae~ 



lim E 

T TO ->0 



Together with ( |A.l 1, this implies E 



[L,L\ 



A 1 



{x 2 A e)F(dx) 



0, because e was arbitrary. As 



T m I 0, [L, L]T m /T m thus converges to zero in probability, and hence a.s. along a 
further subsequence T m ' k ■ This shows that on a set of probability one, the sequence 
[L, L]x m /T m has the cluster point and, moreover, that is the unique cluster 
point. Therefore [L, V\T m /T m — > a.s. and by the right-continuity of [L,L]t the 
claim follows. 

substitute [L, L]r m by L\ . Now note that L is a Levy process 



2.4 



To show Lemma 

with triplet (0, 0, F(dx)) relative to the truncation function h(x) = x, which can be 
used because L is integrable. Again using |T3} Formula (5.8)], the above arguments 
show ^L|, — > a.s. By independence and stationarity of the increments of L, it 
follows that RV£(T) a.s. □ 
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